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IDENTIFICATION  OF  GENE  EXPRESSION  CHANGES  IN  WHOLE  BLOOD  INDICATIVE 
OF  EXPOSURE  TO  CHEMICALS  WITH  DIFFERENT  TARGET  ORGAN  TOXICITY 

Victor  Chan',  Andrea  Stapleton',  Armando  Soto^,  Kyung  Yu^,  and  Nicholas  DelRaso^ 
'Alion  Science  and  Technology,  Corp.,  and  ^Applied  Biotechnology  Branch,  Division  of 
Biosciences  and  Protection,  Human  Effectiveness  Directorate,  Air  Force  Research  Laboratory, 
Wright-Patterson  AFB,  Dayton,  Ohio. 


ABSTRACT 

Toxicogenomic  profiling  is  a  molecular  approach  to  capture  the  global,  transcriptomic  response 
of  a  biological  system  caused  by  perturbations  resulted  from  chemical  exposure.  Coupled  with 
advanced  bioinformatic  techniques,  it  allows  for  the  elucidation  of  the  molecular  mechanisms  of 
chemical  toxicity,  as  well  as  the  identification  of  novel  biomarkers  predietive  for  chemical 
exposure.  In  an  effort  to  investigate  if  gene  expression  changes  in  whole  blood  can  serve  as 
novel  biomarkers  indicative  of  organ-specific  malfunction  resulted  from  chemical  exposure,  we 
have  initiated  a  transcriptomic  profiling  study  using  chemicals  with  different  target  organ 
toxicity  in  a  rat  model.  The  first  three  compounds  chosen  were  a-naphthylisothiocyanate  (ANIT, 
0.1-100  mg/kg),  chlorpyrifos  (CPF,  0.5  -  50  mg/kg)  and  2-bromoethylamine  (BEA,  1.0  -  250 
mg/kg),  which  are  model  hepatotoxin,  neurotoxin  and  nephrotoxin,  respectively.  Gene 
expression  profiles  of  blood  samples  isolated  from  control  and  chemical-treated  rats  were 
determined  using  DNA  microarrays  (Affymetrix  GeneChip,  RAE230A).  The  expression  levels 
of  584,  303  and  441  genes  were  significantly  altered  after  chemical  treatments  as  determined  by 
t-test,  ANOVA  and  SAM  (significance  analysis  of  microarray),  respectively.  Close  examination 
of  these  gene  lists  revealed  that  there  are  a  total  of  1075  unique  genes  that  collectively  comprised 
these  lists.  Of  these  unique  genes,  213  were  identified  by  two  methods,  while  40  genes  were 
common  in  all  three  methods.  Using  a  dataset  of  the  expression  profiles  of  these  40  genes  as 
input  data  in  hierarchical  clustering  and  principal  component  analysis  demonstrated  that  these  are 
useful  discriminators  for  separating  animal  groups  treated  with  different  chemicals.  This  result 
strongly  supports  the  premise  that  whole  blood  is  a  suitable  surrogate  tissue  for  monitoring  gene 
expression  changes  for  detection  of  organ-specific  toxicity.  Future  efforts  will  involve 
transcriptomic  profiling  of  additional  chemicals,  and  validate  gene  expression  changes  identified 
in  this  experimental  model  as  biomarkers  for  specific  exposure  paradigms  and  real-world 
exposure  scenarios. 


OBJECTIVES 

1.  To  investigate  whole  blood  gene  expression  changes  in  rats  exposed  to  a- 
naphthylisothiocyanate,  chlorpyrifos  and  2-bromoethylamine,  the  model  hepatotoxin, 
neurotoxin  and  nephrotoxin,  respectively. 

2.  To  compare  gene  expression  changes  in  rats  exposed  to  chemicals  with  different  target  organ 
toxicity. 


3.  To  investigate  if  whole  blood  gene  expression  changes  in  rats  can  be  used  as  surrogate 
biomarkers  in  rats  exposed  to  chemicals  with  different  target  organ  toxicity. 


METHODS  AND  MATERIALS 

1.  Animals,  Chemical  Treatment  and  RNA  Extraction 

a.  Male  Fischer  F344  rats,  222-258g,  were  obtained  from  Charles  River  Laboratories. 

b.  Animals  were  dosed  by  oral  gavage  with  chemicals  or  vehicle  alone: 

•  ANIT  in  com  oil:  0.1,  0.5,  1.0,  10, 20,  50  &  100  mg/kg 

•  CPF  in  com  oil:  0.5,  1.0, 10,  30  &  50  mg/kg 

•  BEA  in  normal  saline:  1.0,  5.0,  15,  50, 150  &  250  mg/kg 

c.  Control  and  chemical-treated  animals  was  sacrificed  by  CO2  inhalation  at  96  hours  post  dose 

d.  1.5-2. 5  ml  blood  from  each  animal  was  collected  into  a  syringe  and  immediately  added  into 
each  PAXgene  tube.  The  sample  was  immediately  inverted  several  times  to  ensure  complete 
mixing  with  the  reagents.  After  incubation  at  room  temperature  overnight  to  allow  for  lysis 
of  the  blood  cells,  the  samples  were  stored  at  -80°C  until  further  processing.  The  extraction 
of  total  RNA  from  the  processed  blood  samples  was  performed  as  recommended  by  the 
manufacturer. 


2.  Gene  Expression  Profiling  and  Fingerprinting 

a.  Affymetrix  oligonucleotide  GeneChip  arrays  (RAE230A)  were  used  for  mRNA  expression 

(trahscriptomics)  profiling  according  to  the  Affymetrix  GeneChip  Expression  Analysis 
Technical  Manual 

b.  Total  RNA  was  isolated  from  tissues  isolated  from  treated  or  untreated  animals  using  Qiagen 
RNeasy  Mini  kit  (Qiagen,  Valencia,  CA). 

c.  Fifteen  pg  of  total  RNA  was  used  for  first  strand  cDNA  synthesis  using  the  Superscript 
Choice  system  (Invitrogen  Corporation,  Carlsbad,  CA)  in  the  presence  of  an  oligo-(dT)24 
anchored  T7  primer  (Proligo,  Boulder,  CO)  followed  by  2nd  strand  synthesis. 

d.  Biotinylated  RNA  was  synthesized  using  BioArray  high  yield  RNA  transcript  labeling  kit 
(Enzo  Diagnostics,  Inc.,  Farmingdale,  NY)  and  fifteen  pg  of  fragmented  bio-cRNA  was  then 
hybridized  with  Affymetrix  Rat  Expression  230A  gene  array. 

e.  Hybridized  and  stained  chips  were  scanned  with  Affymetrix  GeneChip  Scanner  3000  and  the 
intensity  was  captured  using  Affymetrix  GeneChip  Operating  Software  1.1  according  to 
Affymetrix  standard  procedures. 

£  A  single  expression  level  for  each  gene  was  derived  from  the  11  probe  pairs  representing 
each  gene;  1 1  perfectly  matched  (PM)  and  1 1  mismatched  (MM)  control  probes.  The  MM 
probes  act  as  specificity  controls  that  allow  for  the  direct  subtraction  of  background  cross¬ 
hybridization  signals. 

3.  Data  Filtering,  Statistical  Analysis  and  Visualization  of  Analysis  Results 

a.  The  gene  expression  status  and/or  signal  quality  were  determined  using  the  Detection 
Algorithm  of  the  Affymetrix  Data  Mining  Tool  (DMT  v.  3.1).  Genes  received  “Absent”  call 
across  all  samples  in  the  experiment  of  one  chemical  were  considered  as  “uninformative”. 


Genes  with  uninformative  status  in  two  or  more  experiments  will  be  removed  from  the 
dataset. 

b.  Gene  expression  changes  were  identified  using  three  different  statistical  methods:  t-test, 
ANOVA  and  SAM  using  the  Affymetrix  DMT  and  TIGR  TMeV. 

c.  T-test:  Gene  expression  changes  were  determined  for  each  chemical. 

•  Only  genes  with  p<0.01  in  two  consecutive  doses  are  considered  significant. 

•  All  the  significant  genes  were  then  combined  as  significant  gene  expression  changes. 

d.  ANOVA:  Gene  expression  changes  were  determined  for  each  chemical. 

•  Only  genes  with  p<0.01  are  considered  significant. 

•  p-value  was  determined  based  on  1000  permutations. 

•  False  discovery  rate  for  multiple  measurements  was  controlled  by  adjusted  Bonferroni 
procedure  [8,  9]. 

e.  SAM:  Gene  expression  changes  were  determined  using  multi-class  analysis  procedure. 

•  Four  classes  were  assigned  to  the  samples. 

•  Control:  Control  samples  from  all  experiments  were  pooled. 

•  ANIT:  All  ANIT  samples  assigned  to  this  class  regardless  the  doses. 

•  CPF:  All  CPF  samples  assigned  to  this  class  regardless  the  doses. 

•  EBA:  All  BEA  samples  assigned  to  this  class  regardless  the  doses. 

f.  Classification  of  samples  using  different  datasets  as  input  data  by  Principal  Component 
Analysis  (PCA)  was  performed  using  Partek  Pro. 


RESULTS 

1.  6113,  6299  and  5778  genes  were  identified  as  “informative”  in  the  ANIT,  CPF  and  BEA 
datasets,  respectively.  Of  these  genes,  5426  genes  were  informative  in  all  three  datasets, 
while  599  genes  were  informative  in  two  of  the  three  datasets.  After  the  removal  of  the 
uninformative  genes,  6025  genes  were  used  in  subsequent  statistical  analysis. 

Insert  Figure  1  (Slide  1) 

Figure  1.  Venn  diagram  of  informative  genes  in  each  dataset 


2.  314,  72  and  266  genes  were  identified  as  significantly  changed  in  the  ANIT,  CPF  and  BEA 
datasets,  respectively,  in  t-test.  Of  these  genes,  only  66  genes  were  significantly  changed  in 
at  least  two  datasets,  while  the  majority  of  gene  expression  changes  appear  to  be  chemical 
specific.  A  total  of  584  independent  genes  were  significantly  changed  in  these  three  datasets 

Insert  Figure  2  (Slide  2) 

Figure  2.  Venn  diagram  of  significant  gene  expression  changes  as  identified  by  t-test  after 
chemical  treatments 


3.  92,  94  and  146  genes  were  identified  as  significantly  changed  in  the  ANIT,  CPF  and  BEA 
datasets,  respectively,  in  ANOVA.  Of  these  genes,  only  28  genes  were  significantly  changed 
in  at  least  two  datasets.  This  is  consistent  with  the  results  of  t-test  that  the  majority  of  gene 
expression  changes  are  chemical  specific.  A  total  of  303  independent  genes  were 
significantly  changed  in  these  three  datasets. 

Insert  Figure  3  (Slide  3) 

Figure  3.  Venn  diagram  of  significant  gene  expression  changes  as  identified  by  ANOVA 
after  chemical  treatments 


4.  A  total  of  441  genes  were  identified  as  significantly  changed  multi-class  SAM  analysis. 

5.  Comparison  of  the  gene  lists  obtained  different  statistical  analysis  revealed  that  a  total  of 
1075  independent  genes  were  statistically  significantly  in  the  combined  dataset.  Of  these 
1075  genes,  213  genes  were  identified  in  at  least  two  methods,  while  40  genes  were 
identified  in  all  three  methods. 


Insert  Figure  4  (Slide  4) 

Figure  4.  Venn  diagram  of  significant  gene  expression  changes  as  identified  by  different 
statistical  methods 


6.  Classification  of  all  samples  using  the  expression  profiles  of  these  1075, 213  and  40  genes 
were  performed,  and  the  result  demonstrated  the  expression  profiles  of  the  40  common  genes 
provide  the  best  discrimination  regarding  chemical  treatments,  as  well  as  the  doses  of  the 
treatment. 


Insert  Figure  5  (Slide  5) 

Figure  5.  Principal  Component  Analysis  of  all  samples  using  the  expression  profiles  of  all 
1075  significantly  changed  as  identified  by  different  statistical  methods 


Insert  Figure  6  (Slide  6) 

Figure  6.  Principal  Component  Analysis  of  all  samples  using  the  expression  profiles  of  213 
genes,  which  were  identified  as  significantly  changed  by  at  least  two  statistical  methods 


Insert  Figure  7  (Slide  7) 

Figure  7.  Principal  Component  Analysis  of  all  samples  using  the  expression  profiles  of  40 
genes,  which  were  identified  as  significantly  changed  by  all  three  statistical  methods 


7.  Examination  of  the  functions  of  these  40  genes  revealed  that  major  biological  processes  are 
affected  including  cytoskeleton  organization  and  biogenesis,  ATP  synthesis,  carbohydrate 
metabolism,  cell  adhesion,  protein  biosynthesis,  and  stress  response/xenobiotic  metabolism. 

Table  1.  Description  of  the  40  common  genes  identifled  as  significantly  changed  using  three 
different  statistical  methods 


Insert  Table 


CONCLUSION  and  DISCUSSION 

1.  Exposure  of  rats  to  chemicals  with  different  target  organ  toxicity  results  in  gene  expression 
changes  in  whole  blood. 

2.  Different  gene  sets  were  identified  as  significantly  changed  using  different  statistical 
methods 

3.  Classification  of  samples  grossly  corresponding  to  the  chemical  and  dose  of  the  treatment 
was  achieved  using  the  gene  expression  profiles  of  the  40  genes  that  were  identified  as 
significantly  changed  using  all  three  statistical  methods. 

4.  Although  common  biological  processes  were  affected  by  chemical  treatments  with  different 
target  organ  toxicity,  successful  classification  of  the  samples  derived  from  different 
experiments  suggested  that  the  expression  of  these  genes  were  altered  in  a  chemical- 
dependent  manner. 

5.  The  result  of  this  study  suggests  that  whole  blood  gene  expression  changes  may  be  used  as 
surrogate  biomarkers  for  exposure  to  chemicals,  including  those  with  different  target  organ 
toxicity. 
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Table  1.  Description  of  the  40  common  genes  identified  as  significantly  changed  using  three 
different  statistical  methods 


Probe  Set  ID  |  Gene  Title 


actin  related  protein  2/3  complex,  subunit  4  (predicted) 


1386996_at  myosin  light  chain,  regulatory  B 


1367605_at  profilin  1 


EF  hand  domain  containing  2  (predicted) 


ATPase,  H+  transporting,  VO  subunit  C 


ATPase,  H+  transporting,  VO  subunit  E  isoform  1 


integral  membrane  protein  2C  (predicted) 


glyceraldehyde-3-phosphate  dehydrogenase 


aidehyde  dehydrogenase  2 


similar  to  pyruvate  kinase  (EC  2.7.1.40)  Isozyme  M2 


aldolase  A 


phosphatidylserine-specific  phospholipase  A1 


CD151  antigen 


CD44  antigen 


inositol  polyphosphate-5-phosphatase  A  (predicted) 


fibroblast  growth  factor  (acidic)  intracellular  binding 
rotein 


Nucleosome  assembly  protein  1-like  1 


leukocyte  immunoglobulin-like  receptor,  subfamily  B, 
member  3  (oredicted) 


pro-platelet  basic  protein 


Ornithine  decarboxylase  antizyme  2 


eukaryotic  translation  elongation  factor  2 


eukaryotic  translation  initiation  factor  3  subunit  k 
redicted) 


ribosomal  protein  LI  5 


ribosomal  protein  L7a  (predicted) 


mitochondrial  ribosomal  protein  L41  (predicted) 


peptidylprolyl  isomerase  F  (cyclophilin  F) 


clusterin 


giutathione-S-transferase,  pi  1 


beta-transducin  repeat  containing 


coiled-coil-helix-coiled-coil-helix  domain  containing  3 
redicted) 


similar  to  RIKEN  cDNA  1 1 10059G10  (predicted) 


Transcribed  locus 


Transcribed  locus 


Transcribed  locus 


Transcribed  locus 


Transcribed  locus 


Transcribed  locus 


GO  Biological  Process  Description 


actin  filament  polymerization 


actin  cytoskeleton  organization  and  biogenesis,  protein  import 
into  nucleus.  Rho  orotein  sional  transduction 


cell  motility,  cytoskeleton  organization  and  biogenesis 


Icell  motility,  cytoskeleton  organization  and  biogenesis 


cytoskeleton  organization  and  biogenesis,  muscle  development 


cytoskeleton  organization  and  biogenesis,  actin  polymerization 
and/or  deoolvmerization 


muscle  development 


ATP  synthesis  coupled  proton  transport,  lysosomal  lumen 
acidification 


IATP  synthesis  coupled  proton  transport 


ATP  binding 


glucose  metabolism,  glycolysis,  apoptosis 


carbohydrate  metabolism,  alcohol  metabolism 


Icarbohydrate  metabolism,  alcohol  metabolism 


fructose  metabolism,  glycolysis 


llipld  metabolism 


cell  adhesion 


cell  adhesion,  defense  response,  ureteric  bud  branching 


cell  communication,  inositol  phosphate  dephosphorylation 


signal  transduction 


DNA  replication,  nucleosome  assembly 


immune  response 


platelet  formation,  regulation  of  progression  through  cell  cycle, 
lucose  transDort.  defense  response  to  bacteria 


polyamine  biosynthesis 


I  protein  biosynthesis,  translational  elongation 


protein  biosynthesis,  translational  initiation 


[protein  biosynthesis 


protein  biosynthesis 


mitochondrial  protein  biosynthesis 


[protein  folding 


response  to  stress,  apoptosis 


xenobiotic  metabolism 


[ubiquitin  cycle 


mitochondrial  gene  transcription 


